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Abstract 

Here, we give upper and lower bounds on the count of positive 
integers n < x dividing the nth term of a nondegenerate linearly 
recurrent sequence with simple roots. 
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1 Introduction 



Let {u n } n >o be a linear recurrence sequence of integers satisfying a homoge- 
neous linear recurrence relation 

u n+k = aiu n+k -i H h a,k-iu n+1 + a k u n , for n = 0, 1, . . . , (1) 

where a±, . . . ,a k are integers with a k ^ 0. 

In this paper, we study the set of indices n which divide the corresponding 
term u n ; that is, the set: 

■A4 := {n > 1 : n\ u n }. 

But first, some background on linear recurrence sequences. 

To the recurrence (PQ), we associate its characteristic polynomial 

m 

f u (X) := X k - ai X k - 1 a k ^X - a k = J](X - G Z[X], 

i=l 

where ai,...,a m G C are the distinct roots of f u (X) with multiplicities 
<j\i ■ ■ ■ , cr m , respectively It is then well-known that the general term of the 
recurrence can be expressed as 

m 

Un = ^2Ai(n)a?, for n = 0, 1, . . . , (2) 

i=i 

where Ai(X) are polynomials of degrees at most Oi — 1 for % = 1, . . . , m, with 
coefficients in K := Q[«i, . . . , a m ]. We refer to [B] for this and other known 
facts about linear recurrence sequences. 

For upper bounds on the distribution of J\f u , the case of a linear recurrence 
with multiple roots already can pose problems (but see below). For example, 
the sequence of general term u n = n2 n for all n > having characteristic 
polynomial f u {X) = (X — 2) 2 shows that M u may contain all the positive 
integers. So, we look at the case when f u (X) has only simple roots. In this 
case, the relation (j2J) becomes 

k 

u n = Y, A i<, for rc = 0,l,.... (3) 
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Here, Ai,...,Ak are constants in K. We may assume that none of them 
is zero, since otherwise, a little bit of Galois theory shows that the integer 
sequence {w n }n>o satisfies a linear recurrence of a smaller order. 

We remark in passing that there is no real obstruction in reducing to the 
case of the simple roots. Indeed, let D G N be a common denominator of all 
the coefficients of all the polynomials Aj(Jf) for i = 1, . . . , m. That is, the 
coefficients of each DAi are algebraic integers. Then 

m m 

Du n = DAi(0)a$ + D(Mn) ~ MVX- 
i=i i=i 

If n E Afu, then n \ Du n . Since certainly n divides^ the algebraic integer 

m 

i=i 

it follows that n divides YllLi DAi(0)a^. If this is identically zero (i.e., 
Ai(0) = for all % = 1, . . . ,m), then we are in an instance similar to the 
instance of the sequence of general term u n = n2 n for all n > mentioned 
above. In this case, M u contains at least a positive proportion of all the 
positive integers (namely, all n coprime to D). Otherwise, we may put 

m 

w n = Y DM Q ) a i for n = 0, 1, . . . . 

i=0 

A bit of Galois theory shows that w n is an integer for all n > 0, and the 
sequence {w n } n >o satisfies a linear recurrence relation of order £ := #{1 < 
i < m : Ai(0) ^ 0} with integer coefficients, which furthermore has only 
simple roots. Hence, Af u C Af w , and therefore there is indeed no loss of 
generality when proving upper bounds in dealing only with linear recurrent 
sequences with distinct roots. 
We put 

A u := J] ( ai - a,) 2 = disc(/ u ) (4) 

l<i<j<k 

1 Here, for two algebraic integers a and f3 and a positive integer m we write a = j3 
(mod m) to mean that (a — f3)/m is an algebraic integer. When = we say that m 
divides a. 



3 



for the (nonzero) discriminant of the sequence {w n } n >o, or of the polynomial 
f u (X). It is known that A u is an integer. We also assume that (u n ) n > is 
nondegenerate, which means that on/oLj is not a root of 1 for any 1 < i < 
j < m. Throughout the remainder of this paper, all linear recurrences have 
only simple roots and are nondegenerate. 

When k = 2, u — 0, U\ — 1 and gcd(a!,a 2 ) = 1, the sequence {u n } n > 
is called a Lucas sequence. The formula ([3]) of its general term is 

a n _ a n 

U n = , for 71 = 0,1,.... (5) 

That is, we can take A\ = l/(a.\ — a 2 ) and A 2 = — l/(ai — 0*2) in the 
formula of the general term ([3]). In the case of a Lucas sequence (u n ) n > , 
the fine structure of the set M u has been described in [T5] (see the references 
therein and particularly [8]). We also note that divisibility of terms of a 
linear recurrence sequence by arithmetic functions of their index have been 
studied in [13] (see also [12] for the special case of Fibonacci numbers). 

For a set A and a positive real number x we put A(x) = An[l,x}. 
Throughout the paper, we study upper and lower bounds for the number 
#J\f u (x). In particular, we prove that Af u is of asymptotic density zero. 

Observe first that if k = 1, then u n = Aa™ holds for all n > with 
some integers /I / and a\ {0, ±1}. Its characteristic polynomial is 
fu(X) = X-a x . It is easy to see that in this case #N u {x) = 0((loga;) Lj(|ai|) ), 
where for an integer m > 2 we use u(m) for the number of distinct prime 
factors of m. So, from now on, we assume that k > 2. 

Note next that for the sequence of general term u n = 2 n — 2 for all 
n > having characteristic polynomial f u (X) = (X — 1)(X — 2), Fermat's 
little theorem implies that every prime is in J\f u , so that the Prime Number 
Theorem and estimates for the distribution of pseudoprimes^] show that it 
is possible for the estimate #7V u (x) = (1 + o(l))a;/logx to hold as x — > 00. 
However, we show that $N u {x) cannot have a larger order of magnitude. 

Theorem 1. For each k > 2, there is a positive constant Co(k) depending 
only on k such that if the characteristic polynomial of a nondegenerate lin- 
ear recurrence sequence {u n } n >o of order k has only simple roots, then the 



2 A pseudoprime is a composite number n which divides 2™ — 2. The paper [14] shows 
that there are few odd pseudoprimes compared with primes, while jlOj (unpublished) does 
the same for even pseudoprimes. 
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estimate 

logx 

holds for x sufficiently large. 

In case of a Lucas sequence, we have a better bound. To simplify no- 
tations, for a posititive integer £ we define log^x iteratively as log-L x : = 
maxjlogx, 1} and for t > 1 as \og e x := log 1 (log^_ 1 x). When £ = 1 we omit 
the index but understand that all logarithms are > 1. Let 

L(x) := exp ^A/logxlog 2 a;j . (6) 

Theorem 2. Assume that {u n } n >o is a Lucas sequence. Then the inequality 

holds as x — > oo. 

It follows from a result of Somer [I7[ Theorem 8] that M u is finite if and 
only if A u = 1, and in this case M u = {1}. 

For Lucas sequences with a 2 = ±1, we also have a rather strong lower 
bound on #J\f u (x). Our result depends on the current knowledge of the 
distribution of y-smooth values of p 2 — 1 for primes p, that is values of p 2 — 1 
which do not have prime divisors exceeding y. We use H(x, y) to denote the 
number of primes p < x for which p 2 — 1 is ?/-smooth. Since the numbers 
p 2 — 1 with p prime are likely to behave as "random" integers from the 
point of view of the size of their prime factors, it seems reasonable to expect 
that behavior of U(x,y) resembles the behavior of the counting function for 
smooth integers. We record this in a very relaxed form of the assumption 
that for some fixed real v > 1 we have 

Il(y v ,y)>y v+ °M (g) 

as y — > oo. In fact, a general result from [3, Theorem 1.2] implies that (JHJ) 
holds with any v G [1,4/3). 

Theorem 3. There is a set of integers C such that C C M u for any Lucas 
sequence u with a 2 = ±1, and such that if (jSJ) holds with some v > 1, we 
have 

#K(x) > #£{x) > x* +0 ^ 

as x —J- oo, where 

:= 1 - 1/v. 
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In particular, since as we have already mentioned, any value of v < 4/3 
is admissible, we can take 

?? = 1/4. 

Furthermore, since (jHJ) is expected to hold for any v > 1, it is very likely that 
the bound of Theorem [3] holds with t? = 1 . 

Finally, we record a lower bound on #J\f(x) when a 2 ^ ±1 but A u ^ 1. 

Theorem 4. Le£ {M n } n >o be any Lucas sequence with A u ^ 1. Then there 
exist positive constants C\ and Xq depending on the sequence such that for 
x > x we have 

#N u {x) > exp( Cl (log 2 x) 2 ). 

Throughout the paper, we use x for a large positive real number. We use 
the Landau symbol O and the Vinogradov symbol <C with the usual meaning 
in analytic number theory. The constants implied by them may depend on 
the sequence {u n } n >o, or only on k. We use cq, c%, . . . for positive constants 
which may depend on {u n } n >o. We label such constants increasingly as they 
appear in the paper. 

2 Preliminary results 

As in the proof of [61 Theorem 2.6], put 

D u {x x , ...,x k ) := det(af J )i< ij < fc . 

For a prime number p not dividing a^, let T u (p) be the maximal nonnegative 
integer T with the property that 

p\ \ [ max{l, \N K/Q (D u (0,x 2 , ...,x k ))\}. 

0<x 2 ,-,x k <T 

It is known that such T exists. In the above relation, x 2 , ■ ■ ■ ,x k are integers 
in [1,T], and for an element a of K we use Nk/q{oi.) for the norm of a over 
Q. Since a\, . . . , a k are algebraic integers in K, it follows that the numbers 
N K /q(D u (0, x 2 , ■ ■ ■ , x k )) are integers. 

Observe that T u (p) = if and only if k = 2 and p is a divisor of A u = 
(at - a 2 ) 2 . 
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More can be said in the case when {« n } n > is a Lucas sequence. In this 
case, we have 

\N K/Q (D u (0,x 2 ))\ = |< -aT\ 2 = |A U | 2 K 2 | 2 , x 2 = 1,2,.... 

Thus, if p does not divide the discriminant A u = {a x — a 2 ) 2 = a\ + 4a 2 of 
the sequence {u n } n >o, then T u (p) + 1 is in fact the minimal positive integer 
i such that p | ue- This is sometimes called the index of appearance of p 
in {u n } n >o and is denoted by z u {p). The index of appearance z u (m) can 
be defined for composite integers m in the same way as above, namely as 
the minimal positive integer I such that m \ U£. This exists for all positive 
integers m coprime to a 2 , and has the important property that m \ u n if and 
only if z u (m) | n. For any 7 £ (0, 1), let 

V un = {p : T u (p) < pi}. 

Lemma 1. For x 1 ,y > 2, the estimates 

#{p:T u (p)<y}«-^, #P u , 7 (a;)« 



logy " log-'' 

hold, where the implied constants depend only on the sequence {M n }n>o- 

Proof. It is clear that the second inequality follows immediately from the 
first with y = x 1 , so we prove only the first one. Suppose that T u {p) < y. 
In particular, there exists a choice of integers x 2 , . . . , Xk all in [1, y + 1] such 
that 

p I max{l, \N k /q(D u (0,x 2 , . . .,x k ))\}. 
This argument shows that 

n P I ]I max{l, \N k /q(D u (0,x 2 , ... ,ar fc ))|}. (9) 

Tu{p)<y l<x 2 ,...,w k <y+l 

There are at most (y + l) fc_1 = 0(y k ^ 1 ) possibilities for the (k — l)-tuple 
(x 2 , . . . , Xk). For each one of these {k — l)-tuples, we have that 

\N K/Q (D u (0,x 2 , . . .,x k ))\ = exp(0(y)). 

Hence, the right hand side in (Q is of size exp(0(y k )). Taking logarithms in 
the inequality implied by OH]), we get that 

logp = 0{y k ). 

Tu(p)<y 



If there are a total of n primes involved in this sum and if Pi denotes the ith 
prime, then 

n 

J2^gPi = 0(y k ), 

i=l 

so that in the language of the prime number theorem, 6(p n ) <C y k . It follows 
that p n <C y k and n <C y k /(klogy), which is what we wanted to prove. □ 

The parameter T u (p) is useful to bound the number of solutions n G 
[l,x] of the congruence u n = (modp). For example, the following is [6j 
Theorem 5.11]. 

Lemma 2. There exists a constant 02(h) depending only on k with the follow- 
ing property. Suppose that {u n } n >o is a linearly recurrent sequence of order 
k satisfying recurrence (CQ). Suppose that p is a prime coprime to a^A^. As- 
sume that there exists a positive integer s such that u s is coprime to p. Then 
for any real x > 1 the number of solutions R(x,p) of the congruence 

u n = (mod p) with 1 < n < x 

satisfies the bound 

Ru(x, P ) < C 2 {k) (jj^ + i 

When {w n } n >o is a Lucas sequence, we put 

Q«, 7 = {P '■ z u(p) < P 1 }- 
The remarks preceding Lemma[T]show that #Q„ i7 (x) = #V U)7 (x) + 0(1). 
Hence, Lemma [1] implies the following result. 

Lemma 3. For x > 1, the estimate 



x 2 ^ 



logx 

holds, where the implied constant depends only on the sequence {wn}n>o- 

As usual, we denote by ty(x, y) the number of integers n < x with P(n) < 
y. By [21 Corollary to Theorem 3.1], we have the following well-known result. 

Lemma 4. For x > y > 1, the estimate 

ty(x, y) = xexp(— (1 + o(l))v\ogv) 

uniformly in the range y > (logx) 2 as long as v — )■ oo ; where 

v := (logx)/ (logy). 
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3 The proof of Theorem H 

We assume that x is large. We split the set J\f u (x) into several subsets. Let 
P(n) be the largest prime factor of n and let y := x 1 ^ loglogx . Let 



A/1 (X) 
M 2 {x) 
Af 3 (x) 



{n < x : P(n) < y}; 

{n < x : n<£Ni(x) and P{n) G V u ,i/(k+l)} 
M{x)\ (U^A/K*)) . 



We now bound the cardinalities of each one of the above sets. 
For N\{x), by Lemma HJ we obtain 

#Mi{x) = V(x,y) = xexp(-(l + o(l))vlogv) = o ( — — ) (10) 



logx 



as x — > oo, where 

log X 



log log X. 



logy 

Suppose now that n e ^(x). Then n = pm, where p = P(n) > 
max{y, P(m)}. In particular, p < x/m therefore m < x/y. Since we also 
have p G V Ut i/(k+i){ x / m ), Lemma [1] implies that the number of such primes 

p < x/m is O (^(x/m) k ^ k+1 ^ , where the implied constant depends on the 

sequence {u n } n > . Summing up the above inequality over all possible values 
of m < x/y, we get 



x/y dt 



t k/(k+i) 

l<m<x/y " 1 (11) 

x /y X 



((k + l)x fe/(fe+1) )t 1/(fc+1 



Now let n G M^{x). As previously, we write n = pm, where p = P(n) > y. 
We assume that x (hence, y) is sufficiently large. Thus, m < x/p < x/y. 
Since n G Af u , we have that n \ u n , therefore p \ u n . Furthermore, T u {p) > 
We fix p and count the number of possibilities for m. To this end, 
let {wi}i>o be the sequence defined as we = u v i for all t > 0. This is a 
linearly recurrent sequence of order k. We would like to apply Lemma [2] to 
it to bound the number of solutions to the congruence 

w m = (mod p), where 1 < m < x/p. 
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If the conditions of Lemma [2] are satisfied, then this number denoted by 
Rw{x/p,p) satisfies 

R w (x/p,p) < c 2 {k) - + 1 



pT w (p) 

Let us check the conditions of Lemma [2j Note first that if a%, . . . , are the 
characteristic roots of {« n } n > , then of, . . . , a\ are the characteristic roots 
of {wt}i>i. Hence, 



MX) = ]J(x 



In particular, the term a w ^ corresponding to the recurrence {wi}i>i satisfies 
Q>w,k = assuming that y > 2. Thus, assuming further that y > |ofe|, we 
then have that p does not divide at, therefore p does not divide a w ^ either. 
Next, note that 



Modulo p, we have that 

A w =l II (ou-aj)] =K (modp). 




From the above congruence, we easily get that p | if and only if p | A u . 
Thus, assuming that x is sufficiently large such that y > \A U \, we then have 
that p \ A u , therefore p \ A w either. 

So far, we have checked that p does not divide a w ^A w , which is the first 
assumption in the statement of Lemma [2j 

Let us check the next assumption. 

Note that since p \ A u , the characteristic polynomial f u (X) of {ue}e> has 
only simple roots modulo p. Since p does not divide the last coefficient 
for the recurrence for {u n } n >o either, it follows that this sequence is purely 
periodic modulo p. Let t p be its period modulo p. It is known that t p is 
coprime to p. In fact, t p is a divisor of the number 

lcm [p l - 1 : i = 1,2,..., A;]. 

Choose some n > such that u no ^ 0. Let x be so large such that y > \u no \. 
Since p > y, we have p \ u no . And since gcd(p, t p ) — 1, there exists an integer 
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s with sp = n (mod t p ). Thus, 

w s = u sp = u no (modp). 

In particular, w s is coprime to p. Hence, for x sufficiently large, the second 
assumption from Lemma |5] holds for the sequence {we}e>o- 

Next we show that T u (p) = T w (p). Observe that this number exists (both 
for the sequence {ue}e>o and {we}e>o) because p does not divide a^- Indeed, 
the claimed equality follows easily from the following calculation: 

D w (x 1 , ...,x k ) = det(ap)i<i ij < fe = (det(a^)) p (mod p) 
= D u (xx, . . . ,x k ) p (modp). 

Since n e J\f 3 (x), we have that T w (p) = T u (p) > p 1 /( fc+1 ). 

Lemma [2] now guarantees that the number of choices for m once p is fixed 

is 



R w (x/p,p) < c 2 (k) ( pl+ J (fe+1) + 



To summarize, we have 



A4(«) < £ »(*> + 

■I 1 r' ^ * 



< c 2 (A;) 



< c 2 (/c) ( ir(x) + x 



dt 



£l+l/(fc+l) 



Af 3 (x)<c 2 (k)[ 7 r(x) + 0[ ( j±^)). (12) 



Therefore 



Comparing ([10]). <|XTp and (jl2]). we get that 

#Jv» < c 2 (fc)7r(a;) + ^— + O ( -^—^ ] (13) 

w ~ w w exp((l + o(l))v logv) yyAfc+i) y v 7 

as x — > oo, where the implied constant depends on the recurrence {« n }n>o- 
By our choice of y as x 1//loglogz , the second and third terms on the right side 
of ffT3"j) are both o(n(x)) as x — >• oo, so we have the theorem. 
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4 The proof of Theorem 
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We divide the numbers n G M u (x) into several classes: 

(i) tf x {x) ■= { n e M u {x) : P{n) < L(x) 1/2 }; 

(ii) Af 2 (x) := {n G Af u (x) : P(n) > L(x) 3 }; 

(iii) A7 3 (_) := AT u (x) \ (M(s) UMW). 
It follows from Lemma _ that 

#A/_(_) < ^(x,L(x) 1/2 ) < 



L ^)l+o(l) 

as x — )■ oo. 

For n G A/*« and p | n, we have n = (mod p) and n = (mod -z_(p)). 
For p not dividing the discriminant of the characteristic polynomial for u (and 
so for p sufficiently large), we have z u {p) \ p ± 1, so that gcd(p, 2_(p)) = 1. 
Thus, the conditions n G A/u, p | n, and p sufficiently large jointly force n = 
(mod pz u (p)). Hence, if p is sufficiently large, the number of n G Af u (x) with 
P(n) = p is at most ^(x/pz u (p),p) < x/pz u (p). 

Thus, for large x, 

_ _ rv* ry> _____ T* 

z u (p)<L(x) z u (p)>L(x) 

The first sum on the right has, by LemmaHl at most L(x) 2 terms for x large, 
each term being smaller than x/L(x) 3 , so the sum is bounded by x/L(x). 
The second sum on the right has terms smaller than x/pL(x) and the sum 
of 1/p is of magnitude log log x, so the contribution here is x/L 

(a) i+°(i) as 

x -+ oo. Thus, #A/_(a;) < x/L{x) 1+oi ^ as x -+ oo. 

For any nonnegative integer j, let Ij := [2 J ', 2 J+1 ). For A/3, we cover / : = 
[L(ai) 1//2 , L(x) 3 ) by these dyadic intervals, and we define dj via 2 J = L(x) ai . 
We shall assume the variable j runs over just those integers with Ij not 
disjoint from I. For any integer k, define Vj t k as the set of primes p G 7j with 
z u (p) G Jfc. Note that, by Lemma CD, we have #Vj,k <C 4 fc . We have 



#am«)<5_E E E ^EE E * 



peT jt kneAfu(x) j k P ev jlk VP «W 

P(n)=p 

-E_L_ ___ ^ u ( p )L(a;) 1 /2a J +o(i)^ 
j k peVj.k v 7 v 7 
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as x — > oo, where we have used Lemma H] for the last estimate. For k > j/2, 
we use the estimate 



< 2~ k 



for x large. For k < j/2, we use the estimate 

V 1 < — = 2 fc " J 
^ p. tt (p) 2^2* 

since there are at most order of magnitude 4 fc such primes, as noted before. 
Thus, 

< 2- j/2 = L(x)- a i /2 . 

We conclude that 

#jV 3 (x) < — — - — /Q — 7tt as x -> oo. 

i 

Since the minimum value of t/2 + 1 / (2t) for £ > is 1 occuring at t — 1, we 
conclude that #J\f 3 (x) < x/L(x) 1+ °^ as a; — >• oo. With our prior estimates 
for #N\{x) and #J\f 2 (x), this completes our proof. 

It is possible that using the methods of [5] and [7] a stronger estimate 
can be made. 



5 The proof of Theorem [3 

Since a 2 = ±1, it is easy to see that the sequence u is purely periodic modulo 
any integer m. So, the index of appearance z u {m) defined in Section [2] exists 
for all positive integers m. Further, by examining the explicit formula (jSJ) 
one can see that for any prime power q = p k we have 

Zu(p k ) I z u {p)p k -\ (14) 

In fact this is known in much wider generality. 
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Now, for any real number y > 1 let 



My := 1cm [m : m < y]. 

We say that a positive integer n is Lucas special if it is of the form n = 
2sM y for some y > 3 and for some squarefree positive integer s such that 
gcd(s, My) = 1 and for every prime p \ s we have p 2 — 1 | M y . Let £ denote 
the set of Lucas special numbers. 

We now show that C C Af u for any Lucas sequence u with a 2 = ±1. To 
see this it suffices to show for any n = 2sM y G C and for any prime power 
q \ n, we have ^(g) | n. This is easy for q | s, since then q = p is prime and 
either = p (in the case p \ A u ) or z u {p) \ p ± 1 (otherwise). And since 

p 2 — 1 | My, we have z u (p) | n in either case. 

If q | 2M y , we consider the cases of odd and even q separately. 

• When q is odd, we have q \ M y so q < y. Write q = p k with p prime, 
so that (Ti~4"l) implies z u (q) \ (p — l)p k ~ 1 ,p k or (p + l)p k ~ l . We have 
pk-l <; ^ anc j if 1 < y ; then 2;(g) | M^. The only case not covered is 
p + 1 > y (so p G (y — 1, y]), k = 1, z u (p) = p + 1. Write p + 1 = 2 J m 
where m is odd. Then 2 J | 2M y and m \ 2M y , so p + 1 | 2M y . Thus, in 
all cases, z u (q) \ 2M y so z u (q) \ n. 

• When q = 2 k is a power of 2 with g | 2M y , then since z„(2) G {2,3}, 
we see from ([13]) that either z u (2 k ) \ 2 k or z u (2 k ) \ 3 ■ 2 fc_1 . Since y > 3, 
in either case we have z u (q) \ 2M y . 

We now use the method of Erdos [4J to show that the set C is rather 
large. For this we take 

\ogx 

V '■= q — ^ an d z := y . 

log log x 

We say that q is a proper prime power ii q = £ k for a prime £ and an integer 
k > 2. 

We define V as the set of primes p such that: 

• p g [y + M]; 

• p 2 — 1 is ^/-smooth; 

• — 1 is not divisible by any proper prime power q > y. 
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Note that if q is a proper prime power and q \ p 2 — 1, then q | p ± 1, 
unless g is even, in which case q/2 | p ± 1. Since trivially there are only 
0(t 1,/2 ) proper prime powers q < t, there are only 0(z;i/ _1//2 ) primes p < z 
for which p 2 — 1 is divisible by a proper prime power q > y. Thus, recalling 
the assumption flB]), we obtain 

#P > 1% y)-y + 0(zy- 1 ' 2 ) = 

provided that x — > oo. 

It is also obvious that for any squarefree positive integer s composed out 
of primes p G V, the integer n = 2sM y is Lucas special. 

We now take the set C v (x) of all such Lucas special integers n = 2sM y , 
where s is composed out of 



log x — 2y 
log z 



distinct primes p G V. Since by the prime number theorem the estimate 
My = exp((l + o{l))y) holds as x — > oo, we see that for sufficiently large x 
we have n < x for every n G C v (x). 
For the cardinality of C v (x) we have 

#£,(*) > ( # r P ) > 

Since 

r = (^ 1 + o(l)) i lQ f X and ^ = (logx)^ 1 ^ 1 ) 

log log x r 

as x — > oo, we obtain #£ v (x) > x l ~ l / v+0 ^ as x — >■ oo. Noting that C 
concludes the proof. 



6 The proof of Theorem |4 

Since A u = 0, 1 (mod 4) and A u ^ 0, 1, it follows that |A„| > 1. Let r be 
some prime factor of A u . Then r k G Af u for all A; > (see [HI pages 210 
and 295]). We let k be a large positive integer and look at u r k+4. By Bilu, 
Hanrot and Voutier's primitive divisor theorem (see [1]), u n has a primitive 



15 



prime factor for all n > 31. Recall that a primitive prime factor of u n is 
a prime factor p of u n which does not divide A u u m for any positive integer 
m < n. Such a primitive prime factor p always satisfies p = ±1 (mod n). 
Since there are at most 5 values of k > such that r k < 30 for the same 
integer r > 1, and since u m \ u n if m \ n, we conclude that u r k+± has at 
least r(r k+A ) — 5 = k distinct prime factors p^r, where r(m) is the number 
of divisors of the positive integer m. Let them be p\ < ■ ■ ■ < Pk- Assume 
that |«i| > |of 2 | - For large n, we have that |«i| n / 2 < \u n \ < 2|«i| n (see [6j 
Theorem 2.3]). If . . . , /3& are nonnegative exponents such that 



log(x/r 



A; log ^ 



then r k+i p^ ■ ■ -p^ k < x is in J\f u (see [TTJ Page 210]), so it is counted by 
#J\f u (x). Hence, 



#ac(x) > n 



log(x/r fe+4 ) 
k log pi 



+ 1 > 



log(a;/r fc+4 ) 



> 



( log(x/r 



fc+4\ 



V2r fc+4 log|ai| 

where the last inequality follows from the mean value inequality 

k 



k fi k Y 
Y[\ogp l < I -^2^ogpi j 

i=l \ i=l / 



< 



< 



\0g(\u r k+4\) 



r + log | oil I + log 2 
k 



/2r fe+4 log|ai| x A 
< V * 

for k >2. In the above, we have also used the fact that \u n \ < 2\ai \ n holds 
for all n > 1 with the choice n := r fc+4 . Let C3 := 2 log |«i|. The above lower 
bound is 

^))V^V(i + o 



> 



logo; 
r fe+4 c 3 

logo: 
r k+4 c 3 



k' 2 



log x 
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provided that 

/ C = o(v / loi^), (15) 
as x —7- oo, which is now assumed. So, it suffices to look at 



log a; 

r fc+4 c 3 



exp (A;log(logx/c3) — k(k + 4) logr) . 



Let A := log(logx/c3). In order to maximize the function f(x) := xA — 
x(x + 4) logr, we take its derivative and set it equal to zero to get A — 
2xlogr — 4 logr = 0, therefore x = (A — 4 logr)/(2 logr) = A/(21ogr) — 2. 
Thus, taking k := LA/ (2 logr) — 2j (so that (TT5"]) is satisfied), we get that 
f(k) = f(x) + 0(f(x)) = A 2 /(41ogr) + 0(A), hence 

#A4(x)>e Xp f (10g( ^^ )2 + 0(loglog,; 
\ 4 log r 

/ (log log x) 2 , nn , 

= exp — h C(loglogx / 

y 4 log r 

which implies the desired conclusion with any constant c\ < 1/(4 logr) 



7 Remarks 

We end with a result showing that it is quite possible for #Af u (x) to be large 
under quite mild conditions. Observe that the sequence u n = 2 n — 2 has the 
property that u\ = 0. Here is a more general version of this fact. 

Proposition 1. Let k > 2 and {u n } n >o be a linearly recurrent sequence of 
order k satisfying recurrence (CQ). Assume that there exists a positive integer 
n coprime to such that u no = 0. Then 

#N u (x) > x/logx, 

where the implied constant depends on the sequence {M n }„>o- 

Proof. Since n is coprime to a^, it follows that {u n } n > is purely periodic 
modulo n . Let t no be this period. Now, let 1Z U be the set of primes p = 1 
(mod t no ) such that f u (X) splits into linear factors modulo p. Alternatively, 
1Z U is the set of primes p such that the polynomial f u (X)(X tn o — 1) splits 
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into linear factors modulo p. The set of such primes has a positive density 
by the Chebotarev density theorem. We claim that 

S u C Af u , (16) 

where 

S u ■= {pn ■ p e Tl u and p > n \A u \}. 
The above inclusion implies the desired bound since then 

4MJx) > #K u (x/n )+O(l) > x/hgx. 

So, let us suppose that p > n \A u \ is in TZ U . Then p = 1 (mod t no ), 
therefore p = 1 + At„ for some positive integer A. Thus, pn = + \not no 
and since {«„}„> i is purely periodic with period t„ modulo no, we get that 

u pno = u no+Xno t no = u no = (mod n ). (17) 

Next, observe that since the polynomial f u (X) factors in linear factors mod- 
ulo p, we get that of = a, (mod p) for all % = 1, . . . , k. In particular, 
a? n ° = a" (mod p) for alii = 1, . . . , k. Since the denominators of the co- 
efficient Ai, % = 1, . . . , k, in (j3J) are divisors of A u and p > |A U |, it follows 
that such denominators are invertible modulo p, therefore Aia^ n ° = Aid!^ 
(mod p) for alH = 1, . . . , k. Summing up these congruences for i — 1, . . . , k, 
we get 

k k 

u pno = ^«r° = = u »o = (mod p). (18) 

i=l i=l 

From the congruences (IT7|) and f|T5]) . we get that both p and n divide u pno , 
and since p is coprime to n , we get that prio | u pno . This completes the proof 
of the inclusion f|T6|) and of the proposition. □ 

The condition that no is coprime to at is not always necessary. The con- 
clusion of Propositon [1] may hold without this condition like in the example 
of the sequence of general term 

Un = I0 n - 7 n - 2 • 5 n - 1 for all n > 0, 

for which we can take hq = 2. Observe that k = 4, 

f u {X) = (X - 10) (X - 7){X - 5)(X - 1), 
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and no is not coprime to (Z4 = —350, yet one can check that the divisibility 
relation 2p | u 2p holds for all primes p > 11. We do not give further details. 

Let A4 u (x) be the set of integers n < x with n | u n and n is not of the 
form pn , where p is prime and u no = 0. It may be that in the situation of 
Theorem [U we can get a smaller upper bound for #Ai u (x) than for #J\f u (x). 
We can show this in a special case. 

Proposition 2. Assume that {u n } n >o is a linearly recurrent sequence of 
order k whose characteristic polynomial splits into distinct linear factors 
in 7L\X\. There is a positive constant c^(k) depending on k such that for 
all sufficiently large x (depending on the sequence u), we have #M u (x) < 
x/L(x) c ^ k \ 

Proof. Let y = L(x). We partition Ai u (x) into the following subsets: 



Mi(x) 
M 2 (x) 
M 3 (x) 



{n G M u {x) : P{n) < y}; 

{n G M. u {x) : there is a prime p \ n, p > y, pT u {p) < kx}; 
M u (x)\(M l (x)UM 2 (x)). 



As in the proof of Theorem [2J we see that Lemma H] implies that j^M.\{x) < 
x/L(x) l ' 2+ °^ as x ->• 00. 

As in the proof of Theorem [H 

* MAx)<< £. 

pT u (p)<kx 

We break this last sum according as p G P^i/^+i) and p $ V Ut i/(k+i), respec- 
tively. Lemma [U shows that #7 ? „,i/( fc +i)(0 < t fc/(fc+1) / log*. Thus, 

EX v ~^ X ./' 



and 



Hence, 



pT u (p) ^ p «!/(*+!) 

psP u ,i/(fc+i) pe7 , Uil / (fc+1 ) 



loffo X 

^W^) - 1^ jyyW+l) < y l/(fc+l) • 

P0 > u,l/(,k+l) 



#M 2 {x) < L(x)1/( fc +1)+o(1) as x^oo. 



19 



Suppose now that n G Ais(x). Let p \ n with pT u (p) > kx. Using as 
before the notation t p for the period of u modulo p, as well as the fact that 
T u {p) < kt p and t p \ p — 1 (since f u splits in linear factors over Q[X]), we 
have 

kx < pT u (p) < kpt p < kp 2 , 

so that p > \fx. Thus, n can have at most one prime factor p with pT u (p) > 
kx. So, if n G A4s(x), we may assume that n = mp where p > \fx > m, and 
P{m) < y. Further, we may assume that u m ^ 0. Since p \ u pm and t p \ p—1, 
we have p \ u m . Now the number of prime factors of u m is 0(m). Since the 
number of n G A^3(a;) with such a prime p \ n is 0(x/(pT u (p)) + 1) = 0(1), 
we have 

#M 3 (x) < m ^ y/x^{y/x,y) = as x oo, 

m<y/x 
P(m)<y 

using Lemma HI 

We conclude that the result holds with c 4 (fc) := min{l/5, l/(fc + 2)}, 
say. □ 

Finally, we note that for a given non-constant polynomial g{X) G Z[X] 
one can consider the more general set 

N u ,g := {n > 1 : g(ri) \ u n }. 

We fix some real y < x 1 ^ 2 and note that by the Brun sieve (see [HI Theo- 
rem 2.3]), there are at most 

N 1 « X (^-) (19) 

values of n < x such that g(n) does not have a prime divisor in the interval 
[y,x l l 2 ]. We also note that for a prime p not dividing the content of g, the 
divisibility p \ g(n) puts n in at most degg arithmetic progressions. Thus, 
using Lemma [2] as it was used in the proof of Theorem [TJ the number of other 
n < x with g{n) \ u n can be estimated as 

p|g(n) 
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Using Lemma [T] for any 7 e (0, 1) and the trivial estimate T u (p) 3> logp, we 
derive 

>^ 1 1 / z k ^ z 1 ^' 

Taking 7 to satisfy 

= {z\ogz) 1/{k+1 \ 

we obtain 

I y < iZ - 1 /(fc+l)(log^)-(*+2)/(*+l). 

Summing over dyadic intervals, we now have 



£ ^«^ 1/( ' +1) a°M)- 



^ y ^-(fc+2)/(fc+l) 

Therefore, 



iV 2 « xy- l ' ik+1 \\ogy)-^l^ + x 1 / 2 . (20) 

Taking, for example, y := (logx) fe+1 , we obtain from ( |T9i) and ( 1201) the esti- 
mate 

#A/^(a;) < m + N 2 « x (^p) ■ (21) 

It is certainly an interesting question to bring the bound ( )2T|) to the same 
shape as that of Theorem [TJ However, the method of proof of Theorem Q] 
does not apply due to the possible existence of large prime divisors of g(n). 
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